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INDUCED  RESONANCE  ELECTRON  CYCLOTRON  (IREC) 
QUASI-OPTICAL  MASER 


Introduction 

This  paper  deals  with  a  variation  of  an  electron  cyclotron  interaction 
which  may,  because  of  its  unique  features,  find  applications  as  a  powerful  and 
efficient  short  wavelength  radiation  source.  In  general,  an  electron 
cyclotron  maser  interaction  takes  place  when  a  resonance  exists  between  an 
electron’s  cyclotron  motion  and  the  Doppler-shifted  electromagnetic  wave 
frequency.  One  class  of  electron  cyclotron  interactions  was  proposed 
independently  in  the  late  1950's  by  a  number  of  researchers^"^  and  is  the 
basis  of  the  electron  cyclotron  maser  radiation  source  (gyrotron) 

Typically,  in  gyrotrons,  the  electrons  are  mildly  relativistic  and  the 

* 

interaction  frequency  is  low,  approximately  equal  to  the  electron's  cyclotron 
frequency.  In  the  gyrotron  mechanism,  the  electrons  experience  primarily 
azimuthal  bunching  resulting  in  an  induced  current  density  which  drives  the 
radiation  field.  The  gyrotron  process  is  also  generally  believed  to  be 
responsible  for  radiation  observed  from  the  earth's  magnetic  poles, solar 
atmosphere^®"^®  and  tandem  mirror  devices. Electromagnetic  cyclotron 
interactions  have  also  been  proposed  as  potential  laser  driven  acceleration 
mechanisms^^’^^  as  well  as  a  possible  mechanism  responsible  for  high  energy 
electrons  observed  in  type  III  solar  burst. 

In  this  paper  we  analyze  an  induced  resonance  electron  cyclotron  (IREC) 
quasi-optical  maser.  Here  the  beam  electrons  are  highly  relativistic  and  the 
interaction  frequency  is  much  greater  than  the  electron  cyclotron  frequency. 
The  intrinsic  efficiency  is  found  to  be  insensitive  to  the  electron's 
energy.  In  this  high  frequency  electron  cyclotron  interaction,  the  transverse 
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as  well  as  the  longitudinal  beam  energy  is  converted  into  radiation.  This  is 
unlike  the  usual  gyrotron  interaction  in  which  primarily  transverse  energy  is 
converted  to  radiation.  By  spatially  tapering  the  external  magnetic  field  an 
induced  resonance  condition  can  be  achieved  which  is  shewn  to  substantially 
enhance  the  overall  operating  efficiency.  The  effects  of  finite  electron  beam 
quality,  fractional  energy  and  pitch  angle  spread,  can  be  minimized. 

The  IREC  maser  mechanism  is  similar  to  that  of  the  cyclotron 
autoresonance  maser  There  are,  however,  a  number  of  important 

differences.  In  the  present  mechanism  the  magnetic  field  is  tapered  to 
provide  induced  resonance.  The  radiation  refractive  index  is  chosen  to 
minimize  the  effects  of  finite  beam  quality,  and  the  number  of  electron 
cyclotron  orbits  within  the  interaction  length  is  kept  small,  i.e.,  typically 
less  than  ten.  Furthermore,  the  interaction  takes  place  within  an  open 
resonator. 

We  derive  a  set  of  one-dimensional,  fully  nonlinear  and  relativistic 
orbit  equations  for  particles  in  a  radiation  field  and  spatially  tapered 
static  magnetic  field.  These  equations  are  used  to  study  a  high  frequency 
electron  cyclotron  interaction  within  a  steady  state  oscillator  configuration. 
In  our  oscillator  model  we  assume  a  large  amplitude  temporally  steady  state 
radiation  field  within  a  resonator.  Electrons  are  injected  into  the 
oscillator  along  a  static  magnetic  field  with  a  longitudinal  velocity 
(parallel  to  3^)  much  larger  than  their  gyrational  (transverse)  velocity.  The 
electrons  enter  the  oscillator  uniformly  distributed  in  both  phase  and 
position.  Under  the  appropriate  conditions  the  average  electron  energy  flux 
(power)  will  be  less  at  the  output  of  the  oscillator  than  at  the  input. 
Assuming  an  ideal  resonator,  i.e.,  no  diffraction  or  ohmic  losses,  the  average 
change  in  electron  energy  flux  in  the  steady  state  just  equals  the  output 
radiation  power  of  the  oscillator. 
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Fields  and  Particle  Dynamics 

To  analyze  the  behavior  of  the  IREC  quasi-optical  maser  we  employ  a 
simplified  model  which  consists  of  a  gyrating  electron  beam  streaming  along  a 
nonuniform  static  magnetic  field.  The  electron  beam  is  propagated  through  a 
plane -parallel  plate  resonator  containing  a  large  amplitude  steady  state 
radiation  field  as  shown  in  Fig.  la.  The  vector  potential  associated  with  the 
radiation  field  within  the  resonator  is  a  circularly  polarized  standing  wave 
of  the  form 

A  A 

A(z,t)  -  2A  sin  (kj^z)  (cos  (aij^t)e.^  sin  (wj^t)ey),  (1) 

where  A  is  the  constant  amplitude,  =  2iiil/L  is  the  longitudinal  wave  number, 
2  2  1  /2 

u).  =  c(kj  +  k  ].  is  the  frequency,  k  «  k.  is  the  transverse  wave 

X*  X(  ^  ^ 

number,  2.  is  a  positive  integer  denoting  the  number  of  half  wavelengths  within 
the  resonator  and  L  is  the  resonator  length.  The  external  nonuniform  static 
magnetic  field  is 

A  A  A 

B^(x,y,z)  -  -  (1/2)(xe^  +  ye^ )3B^(z) /3z  +  B^(z)e^,  (2) 

where  B^Cz)  is  the  slowly  varying  longitudinal  component  of  field.  Since  we 
are  interested  in  the  high  frequency  interaction  of  the  gyrating  beam  with  the 
radiation  we  considered  only  the  forward  propagating  component  of  the  field  in 
(1 )  which  is 

A  A 

A^(z,t)  -  A(3in  «j^(z,t)e^  +  cos  <lij^(z,t)ey} ,  (3) 

where  -  k^^^z  ^  dj^t.  To  obtain  the  fully  nonlinear  and  relativistic  particle 
orbits  we  denote  the  transverse  position  and  momentum  of  a  particle  by 
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(4a,b) 


X  =•  X  +  r  3in  0,  y  »  y  -  p  cos  9, 

s  g 

°X  ■  fgx  *  '“J.  •  fjy  *  ’”1  =1"  9. 


(4c,cl) 


respectively.  In  this  representation  (Xg,  yg)  and  (Pg^^,  Pgy)  denote  the 

transverse  coordinates  and  momentum  of  the  particle's  guiding  center,  r  is  the 

Larmor  radius,  Pj_  is  the  magnitude  of  the  transverse  momentum  and  9  is  the 

momentum  space  angle.  We  assume  that  x  ,  y  ,  P  ,  P  r,  P_l  and  i);  -  0  $ 

s  s  sy 

are  slowly  varying  functions  of  z,  i.e.,  change  little  during  a  gyro  period. 

Using  the  Lorentz  force  equation  together  with  the  fields  in  (2)  and  (3)  and 

the  particle  orbit  representation  in  (4),  the  nonlinear  particle  orbit 

equations  for  ?j_,  P  and  can  be  cast  into  the  following  form 
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(y  «  nU  ) 

- U -  cos  4,  >  ^ 


!!o 

3z  ’ 


3z 


ux  ui 

,  an  §  ^  cos  4.  -  ^ 
z  z 


Aa)/c  +  a  ~ 
c 


(Y  ^  nU_) 


o 


!!o 

3z  ’ 


sin  <p. 


(5a) 

(5b) 


(5c) 


where  -  P±/m^c,  Y  »  (1  +  u£  +  is  the  relativistic  mass 

factor,  a  -  |e]A^/m^c  is  the  normalized  radiation  amplitude,  4»  -  0  +  $  is  the 

relative  phase  between  the  radiation  field  and  particle,  n  -  ck/u  »  1  -  (kj_/k)^/2 

is  the  refractive  index  associated  with  the  radiation  field  Aoj  -  (uj(1  -  n6  ) 

z 

-  *1  /Y)/8  is  the  frequency  mismatch  term  and  £2  -  lelB  /m  c  is  the  electron 

vJ  4  O  '  0  0 

cyclotron  frequency.  The  frequency  mismatch  terra  Aw  will  be  shown  to  play  a  central 
role  in  the  cyclotron  interaction.  The  subscript  i!,  denoting  the  mode  number 


has  been  omitted  from  the  frequency  and  wave  number.  Using  (5a)  and  (5b)  the 
spatial  rate  of  change  of  Y  is  given  by 


dz 


a  _  _  coo  V. 
z 


(6) 


Equations  (5a,b,c)  are  the  nonlinear  orbit  equations  for  a  single  particle  in 

the  tapered  external  magnetic  field  given  by  (2)  interacting  with  a  forward 

going  resonator  field  given  by  (3). 

The  interaction  frequency  is  determined  by  the  resonant  frequency 

condition,  ai-v  k-Sl/Y  =  S  Aai,  and  the  eigenf  requency  of  the 
OZ  0  0  oz  o  =>  M  J 

resonator  u  »  ck(l  +  (kj_/k)  /2].  The  simultaneous  solution  of  these  equations 
yield  the  interaction  frequency 


(1) 


(1 


3  )Y^  - 

OZ  oz 


(a  e  AUJ^) 

0  o _ oz  o _ 

1  +  8^  (1  +  3  )Y^  c 

oz  oz  oz 


2y"  a/y. 

oz  O  0 


(7) 
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where  e  -  1-n  »  (k_|_/k)  /2  and  the  last  expression  in  (7)  assumes 

3  *  1,  Iao)  1  <<  /Y  and  lei  <<  1/2Y^  .  Note  that  this  is  similar  to  the 

conventional  free  electron  laser, where  the  interaction  frequency  is  given 
2 

by  Id  -  2Y  ck  where  k,.  is  the  wiggler  wave  number. 
oz  V 

Before  considering  the  effect  of  a  spatially  contoured  external  magnetic 
field  on  the  interaction  we  first  consider  the  case  where  is  constant.  For 
3g  constant  the  orbit  equations  in  (5)  process  two  constants  of  motion  which 
are 


U  -  nY 
z 


(8a) 


U^  ~  2n^Y/id  -  2aUj_  sin  ij/  =  C^. 


(8b) 


With  the  use  of  the  constant  of  motion  it  is  possible  to  estimate  the 
change  in  electron  energy  necessary  to  produce  a  change  in  ij;  equal  to  ir. 
This  energy  change  is  expected  to  yield  an  estimate  for  the  maximum  energy 
change  an  electron  can  undergo  within  the  resonator.  To  this  end  we  expand 
the  Aai  term  in  the  phase  equation  (5c)  about  its  initial  value  Au)^>  writing 
Aoj  -  Au  +  OAu)/9Y)  (Y  -  T  )  where  the  subscript  o  denotes  the  value  at  z  - 

O  0  0 

2 

and  using  (3a)  we  find  that  OAa)/3Y)  =•  (a>(1  -  n  )  -  nAw  )/U  Since 

o  0  oz 

uj  =  (Si  +  U  Aoi  )/(Y  -  nU  ),  (3Ato/3Y)  becomes 

0  oz  o  0  oz  o 

~  n)Aw  /Si^j 

f  J  Atll  1  O  O  OZ  _ O  O  /  Q  \ 

^3Y  ^o  “  Y  (1  -  n8  )  B  Y  • 

o  oz  oz  o 


The  maximum  change  in  electron  energy  at  the  end  of  the  resonator,  z  =  L,  i 


found  by  setting  (3A(i)/3Y)  (Y  -  Y_)L/c  equal  to  *  ir  giving 


where  N  =  «  L/(2irY  v  )  is  the  number  of  electron  cyclotron  orbits  within  the 
CO  o  oz 

resonator.  The  last  inequality  in  (8)  assumes  Iaw^I  is  small  compared  to 
2 

I  Si  (1  -  n  )/Y  (S  -  n)  I .  Equation  (8)  may  be  used  to  obtain  an  estimate  for 
'  o  o  oz  ' 

the  interaction  efficiency.  Although  may  be  large,  say  on  the  order  of  10, 

2 

the  term  (1  -  nB  )/(l  -  n  )  can  be  chosen,  within  certain  limits,  to  be  large 
oz 

resulting  in  high  interaction  efficiencies.  By  spatially  tapering  the  applied 
magnetic  field  the  cyclotron  resonance  condition  can  be  further  extended 
(induced  resonance)  resulting  in  higher  efficiencies.  The  similarity  between 
(10)  and  the  low  gain  estimate  for  efficiency  in  the  free  electron  laser^^  is 
apparent.  In  the  low  gain  PEL  the  intrinsic  (untapered  wiggler)  efficiency  is 


1/2N  where  is  the  number  of  wiggler  periods  within  the  interaction  region. 


Effect  of  Input  Beam  Quality 

Another  important  property  of  the  IREC  maser  is,  if  the  refractive  index 
n  is  appropriately  chosen,  its  insensitivity  to  a  spread  in  the  beam  electrons 
longitudinal  velocity.  Tne  effect  of  an  electron  velooity  spread  is 
manifested  by  a  spread  in  the  frequency  mismatch  terra  Aw.  The  velocity 
spreads  associated  with  the  input  beam  electrons  must  be  kept  small  enough  so 
that  the  frequency  mismatch  spread  does  not  lead  to  a  significant  spread  in 
the  electron  phase  tj).  The  condition  necessary  to  insure  that  the  spread 
in  4)  at  the  end  of  the  resonator  is  much  less  than  ir  is 


I  (^)  * 


(11) 


where  5v  and  6vj_  denote  the  thermal  spread  in  longitudinal  and  transverse 
electron  velocity  about  v^^  respectively.  Using  the  definition  of 
Aw,  the  inequality  in  (11)  becomes 


(Y  q  6  -  nw  -  Aw„)6v 

o  o  oz  0  z 


Y  q  3  6v 
0  o  oj_ 


<< 


TTV  C 
OZ 


(12) 


Since  lAw  i  <<  Y  q  3  or  nw,  the  coefficient  of  the  5v  terra  in  (12)  can  be 
'  o'  o  o  oz  z 

neglected  if  the  index  of  refraction  is 


n  »  Y  q  3  /w. 

0  o  oz 


With  this  choice  of  n  the  inequality  in  (12)  becomes 


53  « 


2N  3  Y 


(13) 


(14) 


It  proves  convenient  to  measure  beam  quality  in  terms  of  AY/Y  and  A0/9  ;  the 

0  0 

fractional  half -width  average  spread  in  the  relativistic  mass  factor  and  the 

velocity  pitch  angle  9=3/3.  Tne  quantities  aY/Y  and  A9/9  are 

0  Oj_  os  0  0 

positive  and  assumed  small.  In  terms  of  AY/Y  and  A9/9  the  average 

0  0 

2  2 

transverse  beam  velocity  spread  is  AS  =  9^(AY/Y^  +  Y^A9/9^)/Y^  which  replaces 
jss^l  in  (14)  giving 

-2  ^oz 

Ae/9  >  Y  AY/Y  «  - ^ — :r.  (15) 

o  o  0  2N  Y^ 

G  Oj|_  O 

If  tne  inequality  in  (15)  is  satisfied,  the  beam's  finite  spread  in  pitch 

angle  and  energy  will  have  little  effect  on  the  interaction.  For  highly 

relativistic  electron  beams  the  inequality  in  (15)  demonstrates  that  a  beam 

pitch  angle  spread  is  more  detrimental  to  the  interaction  than  a  beam  energy 

spread.  This  point  will  be  further  demonstrated  later. 

The  interaction  frequency  is  o)  =  v  k  +  n  /Y  +3  Aoj  which  upon  using 

oz  0  0  oz  0 

the  index  of  refraction  in  (13)  and  neglecting  Au)  gives 


(*)  >  Y  q  (1  -  s'"  ) 

0  0  0j_ 


To  achieve  a  high  frequency  interaction  and  reduce  the  effect  of  an  electron 


pitch  angle  spread,  see  Eq.  (15),  it  is  necessary  to  have  3  <<  1 . 

oj. 


Discussion  and  Numerical  Results 

In  this  paper  we  are  primarily  interested  in  utilizing  the  high  frequency 
interaction  in  an  oscillator  configuration  to  generate  high  power  levels  of 
millimeter,  submillimeter  and  infrared  radiation  efficiently.  To  this  end  the 
particle  equations  in  (5)  are  numerically  integrated  for  an  ensemble  of 
particles  and  the  average  electron  energy  flux  change  is  obtained.  Our 
analysis  is  based  on  the  idealized  model  shown  in  Fig.  1a,  but  perhaps  a  more 
practical  configuration  for  this  mechanism  is  that  shown  in  Fig.  1b.  In  Fig. 
1b  both  the  electron  beam  and  tapered  magnetic  field  are  directed  at  an  angle 
to  the  resonator's  axis.  The  purpose  of  the  tapered  magnetic  field  is  to 
prolong  the  resonance  condition  for  the  average  electron  and  hence  enhance 
efficiency. 

Figure  2  shows  the  efficiency  (average  change  in  electron  beam  energy 

flux  in  passing  through  the  resonator/input  beam  energy  flux)  as  a  function  of 

axial  distance  measured  in  terms  of  electron  cyclotron  orbits  for  a  particular 

set  of  parameters.  The  parameters  for  this  figure  are  ° 

2 

3  Ao)  =  0.035,  n  =  0.9885  and  a  =  lelA/m  c  »  0.025.  Curve  (a)  is 
oz  o  o  '  '  o 

obtained  from  a  small  signal  analysis  of  the  orbit  equations. Curve  (b)  is 
obtained  by  numerically  integrating  the  orbit  equations  for  an  ensemble  of 
electrons  in  a  uniform  magnetic  field.  For  this  curve  the  maximum  intrinsic 
efficiency  is  -  15?  and  is  reached  after  -  5  cyclotron  orbits.  Curve  (c) 
shows  the  efficiency  when  the  external  magnetic  field  is  tapered.  For  this 
curve  the  magnetic  field  is  held  constant  until  the  electrons  phase  bunch.  At 
this  point  induced  resonance  is  achieved  by  linearly  increasing  the  magnetic 
field  by  -  35?  and  the  maximum  enhanced  efficiency  is  -  44?.  This  particular 
spatial  magnetic  field  variation  is  not  optimum,  but  was  chosen  simply  as  an 
illustration.  Comparing  curves  (b)  and  (c)  we  see  that  by  tapering  the 


magnetic  field  the  peak  efficiency  is  increased  by  a  factor  of  three  over  the 

uniform  field  case.  Tapering  the  magnetic  field,  however,  increases  the 

number  of  cyclotron  orbits  needed  to  reach  the  peak  efficiency  =  11). 

Figure  3  shows  the  peak  efficiency  vs.  normalized  radiation  field 

amplitude  for  various  values  of  the  frequency  mismatch  parameter.  In  this 

figure  y  =  4,  8  =*  1 /Y  -  0.25  and  the  number  of  electron  cyclotron  orbits 

O  Oj_  o 

at  maximum  efficiency  and  the  index  of  refraction  are  for  curve  (a) 

N  =»  20,  n  -  0.9950,  for  curve  (b)  N  »  10,  n  *=  0.9925,  for  curve  (c) 
c  c 

N  *■  5,  n  =  0.9885  and  for  curve  (d)  N  =»  3.5,  n  =  0.98^9.  As  the  frequency 
c  c 

mismatch  parameter  3  Aw  /O  is  decreased,  both  the  peak  efficiency  as  well 

02  o  O 

as  the  number  of  electron  cyclotron  orbits  necessary  to  achieve  the  peak 
efficiency  increases.  Although  this  mechanism  is  insensitive  to  a  beam  energy 

spread,  if  the  index  of  refraction  is  n  =  Y  0  8  /u,  Eq.  (13)»  it  can  be 

O  O  02 

sensitive  to  a  pitch  angle  spread,  see  Sq.  (15).  For  the  parameters  chosen  in 
Fig.  3,  the  beam's  pitch  angle  spread  must  satisfy  A6/9^  <<  1/2N^. 

Figure  A  shows  the  effect  on  efficiency  when  the  injected  beam  has  either 
a  spread  in  pitch  angle  or  a  spread  in  energy.  In  this  figure  the  parameters 
are  the  same  as  those  corresponding  to  the  peak  efficiency  (n  •  44^)  in  curve 

(c)  of  Fig.  2.  In  Fig.  4  the  cold  beam  efficiency  is  reduced  by  50%  by  a 

half-width  beam  pitch  angle  spread  of  -  4?  or  by  a  half-width  beam  energy 
spread  of  22?.  In  this  figure  the  number  of  cyclotron  orbits  needed  to  reach 
peaK  efficiency  was  large,  ■  11. 
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Fig.  I  —  Shows  the  high  freguency  electron  cyclotron  maser  configuration  (a)  an  idealized  model 
consisting  of  a  parallel-plate  resonator  and  electron  beam.  The  individual  electrons  stream  along  and 
gyrate  about  a  tapered  magnetic  field,  (b)  here  a  more  realistic  model  consists  of  an  open  resonator  in 
which  both  the  magnetic  field  and  electron  beam  are  directed  at  an  angle  to  the  resonator’s  axis. 


Number  of  Cyclotron  Orbits,  Nc 

Fig.  2  —  Shows  the  efficiency  as  a  function  of  axial  distance  measured  in  terms  of  electron  cyclotron 
orbit.  Curve  (a)  is  obtained  from  a  small  signal  analysis,  curve  (b)  is  for  a  uniform  magnetic  field  and 
curve  (c)  is  for  a  spatially  tapered  magnetic  field. 


Normalized  Radiation  Amplitude,  a  =  ieiA/moC^ 

3  -  Shows  the  peak  efficiency  as  a  function  of  the  radiation  amplitude  for  various  values  of  the  frequency 
mismatch  parameter.  The  normalized  frequency  mismatch  is  shown  for  each  curve. 


^  a  pitch  angle  spread  or  energy  spread.  In  this  figure 
efficiency  is  44%  for  Afl/tfo  “  Ay/yo  =  0  • 
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